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Abstract: 

Transport at a quantum critical point depends sensitively on the relative magnitudes of 
temperature, frequency and electric field. Here we used the gauge/gravity correspondence to 
compute the full temperature and, generally nonlinear, electric field dependence of the electrical 
conductivity for some model critical theories. In the special case of 2+1 dimensions we are also 
able to find the full time-dependent response of the system to an arbitrary time dependent exter- 
nal electric field. The response of the system is instantaneous, implying a frequency independent 
conductivity. We describe a mechanism that rationalizes the instantaneous response. 

1 Introduction 

Quantum critical points continue to be objects of intense attention for reasons ranging from their 
intrinsic elegance to their potential for economically organizing our understanding of the phase 
diagrams of various strongly correlated condensed matter systems pQ. For example, unusual 
power laws in electrical transport have often been attributed to the proximity of the system to 
various quantum critical points, most notably in the case of the cuprates and the heavy fermion 
compounds. 

By definition, quantum critical points lead to singular behavior in various physical quantities 
and a key implication of this is that one cannot understand their transport properties in terms of 
a set of transport coefficients, such as an electrical conductivity. Instead one needs to compute 
the full functional dependence of currents on fields allowing for variation of the frequency uj, 
temperature T and amplitude E (for the case of electrical transport that we consider in this 
paper). This requirement arises as the various limits of small uj, T and E do not commute. 



1 



For example, it was pointed out in [2] that the current j(u),E) established in response to a 
finite frequency and finite amplitude electric field should in general be expected to be a non- 
trivial function of the the dimensionless combination lu/E z ^ 1+z ^ for a quantum critical point 
with dynamic scaling exponent z. This has the implication that a linear response computation 
of the DC transport, which takes the limit E — > before ui — > should be expected to disagree 
with the answer obtained when the limits are reversed. Indeed, in the case studied in [2] this 
was shown explicitly. This discussion parallels a similar competition between the to — > and 
T40 limits first described in (3J. 

Let us now state more precisely the challenge considered in this paper. Quite generally, 
scaling arguments imply that the electrical current at a quantum critical point in d s spatial 
dimensions characterized by a correlation length exponent v and dynamic scaling exponent z 
takes the form: 

T, E) = E d ~^ Y(E/T 1+1 / Z , u/T) . (1) 

Thus the broad aim of a detailed theory of a given critical point is, in addition to computing u and 
z, to compute the function Y and that is what we address in this paper. Let us note that this is a 
non-trivial task. In principle, in the limit E 1+1 / z <C us,T we can resort to linear response theory 
and the Kubo formula. In practice, as noted in [3], this is not always straightforward. However, 
matters are even trickier when we move away from this limit of small E for now the response is 
generically nonlinear in E with a fractional power and one cannot resort to correlation function 
computations at all [HH]. For example, in the limit ui,T <C E scaling predicts the algebraic form 

d s -l + z 

j <xE i+* . (2) 

Two previous computations that have accessed this nonlinear regime (although not the full 
function Y) by means of a transport equation [2] and a clever choice of a dissipative, Gaussian, 
critical point [3]. In this paper we bring newer tools to bear on this, specifically the AdS/CFT 
correspondence. 

For a large class of quantum critical points with a dual gravitational description in terms of 
the AdS/CFT correspondence the finite T non-linear DC (to = 0) conductivity can be readily 
obtained. We will review this construction in the next section. However, AdS/CFT allows 
us to extract much more information about these critical points than just their DC transport. 
Specifically, in the special case of 2+1 dimensional quantum critical points with relativistic 
scaling (z = 1) we are able to construct a full time-dependent solution that describes the response 
to an electric field E(t) with arbitrary time dependence, both at zero and finite temperature, thus 
obtaining the full function Y which is, remarkably, simply a constant. This builds in the physics 
that the resulting current responds to the external field instantaneously and that is simply linear 
in E irrespective of the amplitude. Equivalently, the conductivity j(u), T, E)/E = a(u}/T,u}/^/E) 
is completely independent of the two variables it could potentially depend on. This finding nicely 
extends the result of [5] where it was shown that, in a slightly different theory with gravity dual, 
the conductivity in the absence of electric fields was independent of frequency and thus of the first 



2 



argument. In that work it was argued that this independence was a result of electric/magnetic 
duality in the gravitational description. Our result indicates that this u independence is much 
more general in theories with a gravitational dual. 

Before proceeding to the details some comments on the potential origins and validity of 
nonlinear response are in order. Note first that according to the work energy theorem the total 
energy of the system in a putative transport steady state increases at a constant rate Ej due 
to Joule heating. It was argued in [2] that for small electric field E this Ohmic heating can 
be neglected, as it is of one higher power in the electric field than the current. As E is a 
dimcnsionful quantity what this really means is that Ohmic losses are negligible up to times 
of the order tE = E^ 1 ^ z+1 \ For a time dependent electric field that is ramped up from zero 
to some finite value E over a time period r one would expect that the system should settle to 
a steady state over a time scale set by r and then remain in this steady state up to time tE 
by which point Ohmic heating will have had enough time to build up a non-negligible finite 
temperature. As long as r <C £e there should be a window in which the scaling solution ([2]) 
applies. 

The microscopic picture proposed in [5] as underlying this steady state solution was a com- 
petition between Schwinger pair production (which increases the available densities of charge 
carrier pairs) and scattering losses. Taking these two processes into account one expects the rate 
of change of the current to be given by (focusing on the relativistic case z = 1) 

| = a E^ +1 ^ 2 - bVEj (3) 

The first term is the Schwinger process, the second term is the scattering loss. The E dependence 
of the relaxation rate is once more determined by scaling alone0. The numbers a and b should 
be of order 1 in a generic theory; however, in the examples we will be interested in where the 
field theory has some large intrinsic dimensionless parameters (such as N c , the number of colors 
and A, the 't Hooft coupling in our main example of a supersymmetric gauge theory), a and b 
can of course scale with these. The resulting current in the steady state is j = a/bE ds ^ 2 . Later 
we will make sense of our AdS/CFT results in the light of these ideas. 

Turning to the balance of this paper, Section 2 contains a review of the method of refs [6)17] to 
calculate DC conductivities in a large class of theories with gravitational dual. We will show that 
these systems do indeed exhibit the scaling solutions predicted by [2] ■ In Section 3 we then focus 
on the 2+1 dimensional system and construct the full solution in response to a time dependent 
external electric field. From this we extract the conductivity. We also discuss implications for 
the validity of the simple microscopic model ((3|). In section 4 we discuss the action of S-duality 
in the 2+1 dimensional case and its connection to the results we obtained. We conclude in 
section 5. 

1 For working out the steady state response to a field it does not matter whether we use the field itself to set 
the scale of the scattering rate as we do here or use the current. At zero field the latter is clearly necessary. Later 
we will use J3jl in the time domain but again it will not matter which form we use. 
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2 DC conductivities via AdS/CFT 



2.1 Basic Philosophy 

The basic objects that can be calculated in AdS/CFT are correlation functions of gauge invariant 
operators. Typically one uses this to extract conductivities as a low-frequency limit of a two 
point function (which is just governed by the bulk propagator) and using a Kubo formula with 
all the limitations described in the Introduction. An alternative strategy developed in jHl[7] is 
to turn on E directly as an external source. In this deformed theory one then simply needs 
to calculate a one-point function, (j x ). For this to be possible the system needs to settle to a 
stationary state. For the "probe brane systems" considered in [6l[7], the energy and momentum 
loss rates are of order N, whereas background energy densities are of order N 2 . So for times that 
are not parametrically large in N the background glue acts as an effective heat and momentum 
sink that lets the system come to a stationary state. In this sense the situation is even better 
than in the generic quantum critical system: the loss rates are negligible not just to time scales 
of order tg but actually all the way to NtE- 

To calculate a one-point function all one needs to do is solve for the near boundary behavior 
of all the fields in the bulk. Close to the boundary bulk fields fall off as a power of the radial 
coordinate r (where r — > oo is the boundary) and so these powers can be read off from linearized 
equations of motion. For AdS^+i Maxwell's equations tell us that gauge fields have two modes, 
one constant and one that falls of as r 2 ~ d . The former is non-normalizable. Its coefficient 
corresponds to the source in the field theory. The variation of the on-shell action with respect to 
the source gives the expectation values. Those depend on the subleading coefficient. For a gauge 
field, this can easily be expressed in terms of the field strength. Let us split the coordinates into 
the field theory coordinates and the radial coordinate r. F^ asymptotically approaches a 
constant whose value directly corresponds to in the field theory. So for the DC conductivity 
we are interested in a solution to the bulk equations of motion where F xt asymptotically reaches 
E, the boundary electric field. To turn on a magnetic field in addition one also needs to turn on 
F xy = B. On the other hand, F r/Jt asymptotically falls off as (j^)r 1 ~ d . So for the conductivity 
all we need to find is F xr to extract j x . In the more general setup of [7] one has E and B turned 
on as well as F xr , F rt (which signals a finite density, F rt ~ (jt)r 1 ~ d ) and then also F ry (since 
with both magnetic field and density non-zero one also gets a Hall current). The most general 
system studied in the literature j8] allows in addition for a non-trivial E ■ B. 

For the simplest DC conductivities (no density, no magnetic field) we are looking for a 
solution to the equations of motion of a gauge field in AdS^+i where F xt = E and F xr are the 
only non- vanishing components of the field strength. For a given E, the leading near boundary 
behavior of F xr gives the current and hence the conductivity. But how does one fix F xr in terms 
of El If the field in the bulk is obeying standard Maxwell equations, F xt — E, F xr ~ (j^r 1 ^' 1 
is a solution for any value of E and (j x ). We need one more boundary condition to relate the 
two. As this is supposed to be a result of the non-trivial IR dynamics this boundary condition 
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has to be a regularity condition imposed in the interior of AdS. For Maxwell alone one would 
for example have to demand (j x ) = for any E as otherwise F xr blows up as r — > 0. At finite 
temperature the minimal value of r is the horizon radius and all solutions seem to be fine. Either 
way, this is clearly not the right physics. Maxwell alone does not describe a consistent dual field 
theory. As the field strength grows in the interior one needs to know the full non-linear equation 
obeyed by the gauge field. Identifying the right IR boundary conditions depends on the details 
of the setup. 

For the class of theories that are described by "probe branes" [HUTU] (we'll give examples 
of this type of construction below) the full non-linear action is given by the Dirac-Born-Infeld 
(DBI) action 

S = -N J drdx d ^-dct{g + F) = -M J drdx d ^f~g^J\ + F 2 . (4) 

For the last step we assumed that F A F = 0, otherwise additional terms are present when 
expanding the determinant of g + F in powers of F. J\f is a normalization constant that is 
completely determined by the tension of the brane and the internal geometry; it is proportional 
to N c , the number of colors in the gauge theory whose gravitational dual we are studying. The 
structure of the DBI action gives solutions for the field strength that involve square roots; a 
generic solution that is real near the boundary will turn complex at some radial coordinate. 

The dynamics underlying this feature of the DBI action are related to Schwinger pair produc- 
tion, this time in the bulk. The DBI action is an effective action describing the dynamics of open 
strings ending on the D-brane. In the limit of large string tension all the excited string modes 
can be integrated out and we obtain an effective action for the zero modes of the string, which 
are the gauge fields living on the D-brane as well as the scalar fields describing its transverse 
position. The cndpoints of the string are charged under the gauge field living on the brane. In 
the presence of a small background electric field the tension in the string overcomes the force 
that wants to pull its endpoints apart. As the latter is independent of the applied electric field, 
there exists a critical electric field beyond which strings ending on the brane become unstable. 
This is the point at which the DBI action becomes complex. If we impose such an electric field 
as an initial condition, the system is unstable. 

But even in the presence of a background E field real solutions can be found if one turns on 
F xr (corresponding to a current (j x ) in the dual field theory) in addition to E. It turns out that, 
at least for time- independent E, there exists a unique value of (j x ) for which the solution remains 
real. This is the IR boundary condition that fixes the conductivity. In the time dependent case 
we are going to exhibit an analytic solution that seems to be manifestly free of pathologies. 
Presumably in this case this is also the unique choice. 

As an example, let's review the conductivity of the D3/D7 system corresponding to a 3+1 
dimensional gauge theory with order N 2 gluons coupled to order N c matter fields, bosonic and 
fcrmionic. The metric we are using in the bulk is that of an AdS black hole 

ds 2 = -h(r)dt 2 + dr 2 /h(r) + r 2 dx 2 (5) 



5 



where h(r) = r 2 — r^/r 2 , and the horizon radius rv, sets the temperature as r;,. = irT. Working 
at zero B, zero density the only field we need to turn on is A x = —Et — H(r), so that F xt = E 
and F xr = H'(r). The action in this case becomes 



S = J\f drdtr 3 ^Jl- F 2 t /(hr 2 ) + h/r 2 F 2 r . (6) 



It is easy to see that constant F xt — E solves the t component of the equations of motion 
that follow from this action. The only non-trivial equation comes from the x component of the 
equations of motion. As A x appears in the action only via its derivatives, we can integrate the 
equations of motion and get 



where at this stage (j x ) simply is an undetermined integration constant. Recall that h(r) ~ r 2 
at large r and decreases monotonically as we decrease r with h(r) = at r — r^. So both 
numerator and denominator under the square root change sign at some value of r. The only 
way for H' to stay real is if numerator and denominator change sign at the same value of r. 
The numerator changes sign at a value r* that is set solely by the external field, r 2 h(r*) = E 2 
and hence = E 2 + (7rT) 4 . For small electric field r* = and the corresponding transport 
phenomena are completely dominated by the horizon, which is supposed to roughly encode the 
physics of the longest wavelength fluctuations. As we increase E, r* moves outwards and starts 
sampling the full geometry. So the non-linear transport depends on the details of the geometry. 
The current (j x ) is now completely determined by demanding that the zero of the denominator 
is at the same point r*, that is Af 2 h(r*)r^ = (j x ) 2 and so 



l =^(E 2 + W) 1/4 E. 



where in the last step we plugged in the value of Af corresponding to a D7 on AdS^ x S 3 . In 
the notation introduced in ([lj this implies the identification, 

1/4 
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(9) 



To compare to the answer quoted in [B] one needs to further rescale E and (j x ) by The 
normalization used in [6] is the natural one where a single charge carrier has charge 1. The 
normalization here avoids prefactors of in the DBI action by assigning charge to the 
basic charge carriers. 

This result already shows the crucial points about the limitations of linearized response: the 
linearized response answer can simply be obtained by linearizing (j x ) in E, so a = — | — - T. This 
vanishes in the zero temperature limit. Going first to zero temperature on the other hand gives 
us a — y/~E. Similar expressions can be worked out at finite density and magnetic field 
and they basically all exhibit this feature. One can also obtain the conductivity from the u> — > 
limit of a frequency dependent 2-pt function via a Kubo formula, as was explicitly done for this 
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system in . This calculation is in perfect agreement with the E = answer we obtained above 
(and therefore also fails to capture the finite answer one gets by setting to = T = at finite 
E). In order to get a frequency dependent conductivity from our approach of calculating 1-pt 
functions in the presence of a finite electric field we need to turn on a time-dependent electric 
field. We will do so in the case of a 2+1 dimensional quantum critical point below. 



2.2 Generalizations 



The calculation described above can easily be generalized to other field theories that have a dual 
in terms of gravity coupled to probe branes. By changing the probe brane we can couple matter 
localized on d s + 1 dimensional defects to the same 3+ld conformal field theory living on the 
D3 branes. One can also change the dimensionality of the glue theory, but typically SYM is not 
a conformal theory in spacetime dimensions other than 3+1, so for applications to the physics 
of quantum critical points we will only concern ourselves with conformal defects coupled to 3+1 
dimensional glue. 

For a probe brane in AT = 4 SYM giving rise to a defect with massless matter living on d s 
spatial dimensions, one finds a DC conductivity 

a = a (E 2 + (ttT) 4 )^- 2 ^ 4 (10) 

where <7o is a pure number that can be calculated in terms of the details of the probe (its tension 
and the volume of the internal space it wraps) [B]. Clearly for any d s ^ 2 one fails to find 
the correct non-linear conductivity when working with zero background electric field and then 
taking the limit of small T, whereas a zero T analysis with finite E does indeed give the predicted 
scaling result. 

This calculation can also be generalized to Lifshitz backgrounds, that is metrics that have a 
scale invariance with a dynamical critical exponent z different from 1. The symmetries require 
the metric to take the form [12] 

ds 2 = -r 2z dt 2 + r 2 dx 2 + dr 2 /r 2 . (11) 

One can study the same DBI action on this background to obtain a conductivity via Ohm's 
As in the D3/D7 example we discussed above, reality of the solution imposes a boundary 
condition at a special radial coordinate r* given in terms of the metric components gu and g xx 

by 

g xx \g tt \ = E 2 (12) 

or in other words r\ +z = E. In [6] the conductivity was shown to be given by 

* = *og x r 2)/2 (13) 



2 Some aspects of the DC conductivity in these backgrounds as obtained by a probe brane analysis similar to 
the one performed here were also presented in |13H15| . 
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for any diagonal metric, where the metric component once more has to be evaluated at r*. This 
reproduces 

a = (TqE 1 ^ (14) 
as predicted by the scaling analysis of [2J. 



2.3 Conductivity of the D3/D5 system 

As we have seen above, the pattern of vanishing conductivities when working with E = 0, T — > 
missing the correct zero temperature answer is generic with one exception: 2+1 dimensional 
defects. In that case there is no non-linear dependence of the conductivity on the electric field 
and the limits of zero T and zero E commute. In fact, the conductivity is completely independent 
of both. The easiest example of a 2+1 dimensional defect coupled to the D3 brane glue is given 
by a D5 brane probe on AdS4 x S 2 , as this theory preserves half of the supersymmetries of the 
glue and so stability is guaranteed. In this case (at zero density, zero magnetic field and for 
massless flavors) the conductivity a is given by 

2N f N c 



(15) 



This formula is true at all values of E and T and is independent of the dimensionless ratio E/T 2 . 
Scale invariance alone would allow a non-trivial dependence on this ratio, so the independence 
seems to be a special property of field theories with a probe brane gravitational dual. 

A very similar phenomenon in theories with holographic dual has been observed before in [5] . 
There it was shown that in 2+1 dimensional theories with gravitational dual in linearized re- 
sponse (that is at E — 0) a, which should be a non-trivial function of w/T, is in fact completely 
frequency independent. Together with the result on the E/T 2 independence of the DC con- 
ductivity this strongly suggest that a(uj 7 T,E) is actually constant. In order explore this in our 
framework we want to calculate the one-point function of j x not in response to a constant elec- 
tric field, but to a time varying electric field. We will address this problem in the next section. 
Before we embark on that calculation, let us note that with cither a finite density of charge 
carriers or a finite magnetic field turned on a is no longer independent of E and T and the 
linearized response once more fails to yield the desired zero temperature answer. 

3 Time dependent solutions 

3.1 General Strategy 

For a time-dependent electric field one can still define a conductivity via Ohm's law 

(j x (t)} = [ dra{r)E{t-r). (16) 
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In order to calculate this time dependent (j x (t)) we have to redo the calculation of the world- 
volume gauge fields in the presence of a time dependent electric field. That is the boundary 
condition in the UV of the AdS geometry is that F xt = E(t) at large r. The 2nd condition 
in order to uniquely specify the field strength is still expected to be the requirement that the 
solution is real. 

3.2 Exact time dependent solution in 2+1 dimensions 

So far we have been using the metric in its standard black hole form 

ds 2 = -h(r)dt 2 + ^ + r 2 dx 2 . (17) 
To study time dependent solutions a better suited time coordinate is the incoming v 

t = v + f(r) (18) 

where 

f(r) = -l/h(r). (19) 

At large r we have t = v and both correspond to the field theory time. Causality demands that 
information about the time varying electric field can at best propagate inwards at the speed of 
light. Lines of constant v are spacetime points that are connected by lightrays to a boundary 
event at t = v and so it is natural to look for solutions in the bulk whose time dependence is 
just governed by v, In terms of v the metric reads 

ds 2 = Idvdr - hdv 2 + r 2 dx 2 . (20) 

In this coordinate system the non- vanishing components of the inverse metric are g rv — g vr = 1, 
g rr = h and g xx = 1/r 2 . 

Denoting quantities in the t, r coordinates with bars, we can work out the effect of the change 
of coordinates on field strength in order to connect to the dictionary we established relating F xt 
and F xr to E and (j x ) respectively The solutions we are interested in have only F xv — E(v, r) 
turned on. So F xr — E/r 2 , hence F xr — E/r 2 , F xt = —E/(hr 2 ) and last but not least 

F xr = E/h, F xt = E. (21) 

So F xv — E(v,r) alone is sufficient to encode both a background electric field and a current pro- 
portional to E. Alternatively, one can directly calculate the conductivity in the v, r coordinate 
system from variation of the action with respect to F xr . 

3 Explicitly, the change of variables gives 

pxv _ pxt _|_ pxr pxr _ pxr 

and hence 

Fxr = Fxr — Fxt/h, F xv = Fxt- 
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The Lagrangian for a D5 brane with these worldvolume fields turned on is 



£ ~ V- detfo + T) = r 2 C = rVr 2 (2F lr F ra + hF*~) + 1. (22) 



The resulting equation of motion is 



<>,[ — + - hFxr )+d v r-f) =(. (28) 



together with the Bianchi identity 

We can find one simple analytic solution for any external forcing function E(v) and temperature 
T: 

F xr = 0, F xv =E{v), C = l. (25) 

In the field theory this corresponds to a state with j(t) — aE(t). The solution is real and 
completely well behaved and so presumably corresponds to the unique physical solution. 

What we have thus established is that in 2+ld one can again find (j x (t)) analytically even 
in the presence of a time dependent electric field with the result that 

(j B (t)) = aE(t) (26) 

that is the current adjusts itself instantaneously to the time varying electric field. In terms of the 
general cr(r) this means that cr(r) ~ S(t) and so its Fourier transform <j(uj) is ^-independent. 
So indeed <t(uj, T, E) in the case of a 2+1 dimensional quantum critical point with a gravity plus 
probe brane dual is completely independent of its arguments confirming our earlier suspicion 
based on the independence of E/T 2 at ui = and oiuj/T and E = 0. 

Last but not least, one can calculate the energy density deposited into the field theory using 
the techniques of [TB] (that is one simply calculates the stress tensor of the probe D5) and finds 
that the energy density e is given by 



dtj{t)E{t) (27) 



as it should be. 



3.3 A microscopic mechanism 

Having worked out the conductivity in the holographic theory and having established that, for 
d s = 2, the conductivity is independent of T, E and w, we can go back and wonder whether 
this feature can be understood, at zero temperature, in terms of the microscopic mechanism of 
Schwinger pair production competing with scattering losses. For the case of d s = 2 our model 
equation reads 

f = aE 3 / 2 - b^Ej. (28) 
at 
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From our knowledge of the conductivity in the stationary state we find again the ratio of a and 

*fl,-a-™&. (29) 

7TV A 

To pin down the scaling of a and b with N c and A we need more information. Note that it is 
natural for the Schwinger process to scale with the number NfN c of quarks that can run in the 
loop, so presumably in terms of the N c counting a ~ N c and hence b ~ 1. But we also need to 
understand how a and b scale with the other large parameter, the 't Hooft coupling A. For what 
follows we will assume that b itself is large. A plausible scaling is to take b ~ vA and hence 
a ~ 1 as far as the scaling with A is concerned. This would imply that unsupported currents 
decay with a time scale that scales as tg/vA- But as we will demonstrate momentarily, any 
parametrically large b is consistent with the frequency independent conductivity we found. 

We can solve the 1st order transport equation (|28|) for general time dependent E: 

j(t) = a[ dt' £(t')3/2 e -«r(t,t') (go) 

J — oo 

where the kernel K(t,t') is given by 



K(t,t') = J sjE{t)dl (31) 

For large b the integral in (|30[) is dominated by the upper integration boundary as all earlier 
times are exponentially suppressed. As always when an integral is dominated by an endpoint, 
we can evaluate the leading contribution from integrating by parts 



t 

3K) -b(drK(t,t')f [t) 



0{b~ 2 ) 



= - b E{t)+0(b- 2 ) = aE{t)+0{^) , (32) 

which then makes sense of the instantaneous response. 

It is interesting to note that the considerations embodied in (|28|) also imply that absent an 
electric field the current must relax to zero as 1/t 2 . [The reader should revisit Footnote 1 if this 
is not clear]. Clearly, our result thus far implies that the current instantaneously drops to zero 
once the field is switched off which is a trivial example of this scaling. It is, however, interesting 
to search for but fail to find a constant current carrying solution to get at this result a different 
way that is consonant with the AdS/CFT ideas we have used in this paper. In fact a potential 
such solution, for the case of zero temperature, has already been constructed in [17] . That work 
was mostly concerned with the phase diagram at finite density, but it was noted that in the end 
the answers only depend on {j^j^) as it should be by boost invariance and so one can easily 
extend the calculation in [17] to the case of a spatial current density. For a solution with a 
persistent current in the x direction the bulk solution reads 

p _ O's) (33^ 
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When (j 2 ) = —d 2 + (j'l) is negative, this is a perfectly regular solution. However for the case 
of interest with d = and constant (j x ) there is a singularity at a finite value of r. Clearly the 
persistent current is not a stable background. 

4 S-duality 

Above we have shown by explicit construction of a time-dependent solution that in 2+1 dimen- 
sional field theories with a probe brane dual the conductivity does not depend on either the ratio 
U)/T or us 2 /E. This result extends the earlier findings of [5], who showed that for the holographic 
dual of the M2 brane theory in the absence of an electric field the conductivity is independent of 
cj/T. As this result was obtained via a standard Kubo formula, only the quadratic Maxwell part 
of the bulk gauge field action was important. It was argued in that work that the underlying 
reason for the u) independence was S-duality. 

We found a similar frequency independence for a system governed by the DBI action. Re- 
markably, the DBI action is known to also be invariant under a non-linear generalization of 
S-duality [l"8rJ2"0] . In terms of E, B and D = 5C/5E and H = -SC/5B its action is just the 
standard electro-magnetic duality transformation of Maxwell's equations in matter. In this sec- 
tion we want to see to what extent S-duality is related to the frequency independence in the 
system we analyzed. 

It was argued in [21] that for any 2+1 dimensional conformal theory with a global U(l) 
symmetry one can define a T-transformation and an S-transformation that takes the theory 
into a different conformal field theory (CFT) whose conductivities are related to the ones in the 
original CFT. For CFTs with a bulk dual these operations correspond to an electric-magnetic 
duality transformation. The S-operator acts on the boundary field theory quantities as 

jx -> -<ToE y , j y -> -a E x , E y -> j x /<Jo, E x -> j v /a , B -> -d/a , d -> a Q B (34) 

These transformation rules can be derived most easily from the bulk point of view. The near 
boundary behavior and hence the transformation of the field theory quantities is simply given by 
the standard Maxwell duality relations in AdS4 . Near the boundary all components of the field 
strength approach a constant whereas the components of the inverse metric go to zero, so the 
linearized equations of motion derived from the Maxwell action hold. The above transformation 
rules then follow recalling that with both legs along the field theory directions approaches 
the corresponding field theory quantities, whereas the mixed components^ F Ufl approach ctoJ/j- 

The S-operator can be applied irrespective of whether S-duality is a symmetry of the bulk 
action. In our case S-duality is indeed a symmetry of the classical equations of motion of the bulk 
theory so the theory generated by the S-transformation will have identical transport properties^. 

4 For the purpose of discussing S-duality it is convenient to use a radial coordinate u = 1/r. 
°Thc application of S-duality in the bulk is complicated by the fact that S-duality does not commute with 
the g s — > 0, 4irg s Nc = A fixed 't Hooft limit. S-duality takes the probe D5 brane into a probe NS5 brane. While 
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Let us demonstrate this in detail. According to the results of [7] the conductivities in the D3/D5 
system in the presence of both electric and magnetic background fields as well as with a finite 
density of charge carriers are 



3x 



9xx + B 2 

dB , . 

* = ^+W E (35) 



y/cr*(g* x + B*) + <PE 



where 



9lx = \ (W) 4 + E 2 -B 2 + yj - B*y + 2(^)4(^2 + B 2 ) ~ (7rT)8 j . (36) 

While we already know the transformation properties of the field theory quantities j^, B and 
E, the transformation rules for g xx need to still be determined. g xx stands for the value of the 
spatial components of the metric evaluated and a special radial location, the worldsheet horizon 
it*. As the latter depends on the full u dependent profile of the bulk field strength and not just 
the near boundary behavior, one needs to take into account the S-duality of the full non-linear 
DBI action. From the expressions derived in [7] it can easily be seen that g xx is invariant under 
the S-transformation. The location of the worldsheet horizon it* is determined in [7] as the 
simultaneous vanishing locus of 3 functions £, x ancl a - & X an d a are an expressed in terms 
of z- independent boundary data (currents, electric fields, magnetic field, density) and so their 
transformation under S-duality is determined by the transformation properties of the boundary 
data It is easy to see that under S-duality 

£ -> x/°~o, X -> cto£, a. (37) 

As z* is defined as the vanishing locus of £, \ and a it remains invariant and so does g X x{z*)- 



Applying the transformation (|34|) to the expressions for the conductivities ([35)) one sees that 
in the dual field theory only j y is turned on, whereas the electric field now has both E x and E y 
components. The interesting physical quantities are the magnitude j of the induced current as 
a function of the magnitude of the electric field, as well as the Hall angle, the angle between 
j and E. After some algebra one finds that in the dual theory these two quantities are in 
fact identical to the original ones; electric field and resulting current in the S-dual system are 
completely equivalent to the original one up to an overall rotation in the x-y plane. So the 
S-transformation generated a new theory which, by virtue of the S-duality of the DBI action, 
has identical transport properties. 

However, unlike in the case of the M2 brane gauge theory studied in [5], S-duality does 
not impose any restriction on the conductivity of the original theory. The ui independence of 



both are described by a DBI action, the tension, which appears as a prefactor of the action, is of course different. 
As the classical equations of motion are independent of the prefactor of the action, the classical gravity + probe 
system has equations of motion that are invariant under S-duality and a subsequent constant shift in the dilaton 
which takes the theory back to weak coupling. So we expect the expressions for the conductivities to actually be 
invariant under S-duality. 
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the conductivity therefore seems to not be intrinsically connected to S-duality. In fact, had 
we taken any Lagrangian C(F 2 ) being an arbitrary functioij^] of F 2 in the same AdS or AdS 
black hole background, we would still get that the time dependent (|25|) solves the equations of 
motion for any time-dependent electric field on the boundary, with an instantaneous response of 
the current to the field. However S-duality would no longer be a symmetry of the bulk action. 
The instantaneous response therefore is a genuine property of any field theory with a classical 
holographic dual of the form spelled out in footnote [6] and does not seem to be deeply connected 
to S-duality. 



5 Conclusions 

We have analyzed electrical transport at the strongly coupled quantum critical points which 
have a dual holographic description in terms of probe branes. We find a stationary solution at 
an arbitrary temperature with a constant current and a generally non-linear dependence on the 
electric field. This adds to the small [2|4| catalog of quantum critical points where the non-linear 
dependence has been explicitly computed. Of the previous computations, this present work is in 
the spirit of [2] where a regime of non-linear response is explicitly exhibited for a non-dissipative 
field theory in a large N limit. Indeed, the explicit mechanism involved in |2| provides a heuristic 
guide to our results. The non-linear dependence on E is uniquely governed by scaling properties; 
the coefficient is finite and calculable in our models. Clearly this phenomenon is beyond the 
reach of linear response and so the application of Kubo formulas to the problem of transport at 
critical points misses this effect. 

For the special case of 2+1 dimensions we were able to also find an analytic solution in the 
case of a time-dependent finite electric field, both at zero and finite temperature thus providing a 
fairly general computation of the critical electrical transport. The resulting form is quite unusual 
in that the current instantaneously adjusts to changes in the external field and is always linear in 
the electric field. This extends earlier findings regarding the independence of the conductivities 
in 2+1 dimensional holographic models from ui/T at E — and from E/T 2 at cj = to the full 
lo/T, lu 2 / E plane. While it was earlier argued that this frequency independence is a consequence 
of S-duality, in our case there seems to be no causal connection. While our system does exhibit 
a nice action of S-duality, this does not seem to force an u) independent conductivity all by 

6 By a Lagrangian of the form C(F 2 ) we mean a Lagrangian that is v / — g times an arbitrary function f(F 2 ) 
of F 2 = FnvgW g" v F^i v i . The DBI itself is manifestly of that form with / = Vl + F 2 as can be seen from eq. 
(0). In completely parallel with the DBI equation of motion from eq. l|23|l the equation of motion in the general 
case is 

9r [f'(F 2 )(F xv + hF xr )] + d v [f'{F 2 )F xr ] = (38) 

which is clearly still solved by F xr = (and hence F 2 = 0) with r independent F xv = E(v). Of course, generic 
higher derivative interactions as they would be produced by quantum corrections would change the conclusions. 
In particular couplings of the form RF 2 in the action would spoil the solution J25}- Also, if a background scalar 
field <f> is turned on with a non-trivial profile, couplings of the form F 2 times a function of <f> would ruin the time 
independent solution. 
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itself. Furthermore, we are able to construct a large class of different 3+ld bulk theories that 
all exhibit frequency independent conductivities without being S-duality invariant. 

Altogether, the results in this paper add to the growing set of results that illustrate the 
power of the AdS/CFT correspondence to illuminate issues of principle in the study of quantum 
critical phenomena. 
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